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We present a theoretical study of a hybrid circuit-QED system composed of two semiconducting 
charge-qubits confined in a microwave resonator. The qubits are defined in terms of the charge states 
of two spatially separated double quantum dots (DQDs) which are coupled to the same photon mode 
in the microwave resonator. We analyze a transport setup where each DQD is attached to electronic 
reservoirs and biased out-of-equilibrium by a large voltage, and study how electron transport across 
each DQD is modified by the coupling to the common resonator. In particular, we show that the 
inelastic current through each DQD reflects an indirect qubit-qubit interaction mediated by off- 
resonant photons in the microwave resonator. As a result of this interaction, both charge qubits 
stay entangled in the steady (dissipative) state. Finite shot noise cross-correlations between currents 
across distant DQDs are another manifestation of this nontrivial steady-state entanglement. 

PACS numbers: 



I. INTRODUCTION 

Recent technological progress has made it possible to 
coherently couple superconducting qubits to microwave 
photons on a superconducting chipP This so-called cir- 
cuit quantum electrodynamics (circuit-QEDpl has paved 
the way to new research directions beyond standard cav- 
ity QED systemsP^l Apart from the high degree of tun- 
ability in circuit-QED, most of the novelty comes from 
the fact that the coupling between qubits and microwave 
photons can reach values well above the ones between 
natural atoms and photons in optical cavities. 6 

An interesting alternative to the above ideas is to im- 
plement hybrid circuit-QED 7 with qubits defined in semi- 
conducting quantum dots (QDs)p"^ as has been ex- 
perimentally implemented recent lyPSHIZ! In these hybrid 
structures, the semiconducting QDs are typically coupled 
to normal electronic reservoirs such that electronic trans- 
port may be used to characterize/modify the properties 
of the circuit-QED system. This possibility has remained 
largely unexplored, except for some works analyzing the 
transport-induced lasing states in the resonatorPSED 

In this context, we here analyze how the coupling to a 
common photon mode generates entanglement between 
distant charge qubits realized in double quantum dots 
(DQDs) and how this entanglement manifests in the 
transport properties of the system. In particular, we 
present a detailed analysis of how the electron currents 
across each DQD are modified due to the interaction with 
the photons in the circuit. The coupling of each DQD 
to a common microwave resonator generates an indirect 
coupling between DQDs which gives rise to positive shot 
noise cross-correlations between distant currents across 
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FIG. 1: (Color online) Schematics of the two charge qubits 
coupled to a transmission line resonator. An excess charge in 
each double dot (formed by the Li and the FU dots) defines 
the states of the qubit. Both qubits are attached to electronic 
reservoirs, via the rates Tl,i and Tr.^, such that an electrical 
current pass through them. The qubits are located at the 
ends of the resonator in order to enhance the coupling with 
the electromagnetic field. 



them. We analyze this physics in terms of an effective 
model and show that off-resonant photons are responsible 
for the induced indirect coupling. Moreover, we demon- 
strate that both charge qubits are entangled in the steady 
(dissipative) state due to this cavity-induced coupling. 

The paper is organized as follows. In section [TT] we de- 
scribe the model for two double quantum dots coupled 
to a microwave resonator as well as the master equation 
that governs the dynamics of this open quantum system. 
In section |III| we discuss the stationary transport prop- 



erties (mean value of the current and shot noise) of the 
system. This section is divided in two parts. The first 



part (subsection III A ) reviews the case of a single double 
quantum dot. We then turn to the analysis of the two 
double quantum dot system (subsection III B ) by calcu- 
lating also shot-noise cross-correlations between distant 
currents across each double quantum dot. In Section [TV] 
we focus on the generation of qubit-qubit entanglement 
induced by the common coupling to a microwave pho- 
ton mode, and compare it with the results obtained for 
the cross-correlations in the previous section. In particu- 
lar, we present a detailed analysis of steady-state Bell 
state occupations and demonstrate that indeed cross- 
correlations between distant currents constitute an indi- 
cator of non-local qubit-qubit entanglement. In sectionfV] 
we extend our study to the case of asymmetric couplings 
between each double quantum dot and the microwave 
resonator. Our conclusions are presented in section |VI| 



II. MODEL 

We consider the coupling of the charge states of two 
uncoupled semiconductor DQDs to an electromagnetic 
resonator with a high Q-factor, as for instance the su- 
perconducting transmission line described in the recent 
experiments of Ref. [T5] We assume that the DQDs are 
placed at the ends of the resonator, as schematically de- 
picted in figure [l] In the following, we consider that the 
charging energy on each DQD is the largest energy scale 
of the problem such that, for each individual DQD, an 
excess electron defines the two states of a charge qubit, 
\Li) and \Ri) (i = 1,2), see e.g. Ref. 1221 In this basis, 
the Hamiltonian describing the DQDs reads 



Hpi = 



Z^ 



1 i 



+ tit7„ 



(1) 



where the energy detuning in each DQD is given by et, U 
is the tunneling coupling between dots of the i-th DQD 
and <jj is the j-th Pauli matrix acting on the charge basis 
of each qubit, namely <j\ = \L.i)(Li\ — \Ri)(Ri\ and a % x = 

The transmission line resonator is modeled as a quan- 
tum harmonic oscillator H res — huj r a^a, where a^ (a) 
is the creation (annihilation) operator of photons in the 
resonator with frequency ui r . The charge states of each 
qubit are coupled to the same mode of the resonator, 
such that the coupling term reads 



H,, 



~Y^hg i a-l{a t 



(2) 



Experimentally, typical photon frequency values range 
from 1 to 10 GHz, whereas couplings strengths g ~10- 
30MHz have been rep orted for a single DQD coupled to 



Finally, we consider that each DQD (i=l,2) is coupled 
to electronic reservoirs described by the Hamiltonian 



H, 






leads 



r J 



2^ ^ S kA L ,i C ^,i + e k,i C k R ,i C kH,i} ( 3 ) 



in which c k ^(cfc-j) is the creation (annihilation) opera- 
tor of electrons in the left/right contact, (3 £ L,R, with 
energy et .. The coupling of each DQD to the leads reads: 

H int = ]T E{ y M c U,A>* + h - c - + L ^R) (4) 



where d^mi (dt , R A creates (annihilates) an electron in 

the left/right QD of each DQD, and V ki are the tunneling 
matrix elements. Due to this coupling to the reservoirs, 
situations in which either of the two DQDs (or both) are 
empty need to be considered and hence the Hilbert space 
in the charge sector is spanned by the states |c*i, 0:2)1 
with a = L,R, 0. This transport model can be easily 
extended to a system consisting of several qubits, see 
e.g. Ref. 1231 and is the single-mode version of previous 
studies focusing on bath-mediated interactions! 2 ^- 

The total Hamiltonian of the system is given by H to t — 
H e i + H res + H e - res + Hi ea d s + H int . The dynamics of 
the resonator and the DQDs is described by the master 
equation for the reduced density matrix pit) obtained 
after tracing out the reservoirs degrees of freedom and 
applying a Born-Mar kov a pproximation with respect to 
the Hamiltonian H int 125126] j n ^ e Schrddinger picture the 
master equation reads pit) — Cp with the Liouvillian: 



Cp = -i [H el + H res + Fe-res, p(t)] 
i 



-- (ap{t)a) - 2a t ap(t) + p{t)a*a) 



(5) 
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with the tunneling rates to reservoirs rf = 

^ n J2k i iy& i\ 2 3( £ i,p — £ k,i,p) and where we considered 
the limit of infinite source-drain voltage, pl ~> oo and 
jitjj — > — oo (such that the Fermi functions in the reser- 
voirs become /l = 1 and fn = 0). In this limit, the Born- 
Markov approximation with respect to the coupling to 
reservoirs is essentially exact and, more importantly, the 
physics no longe r dep ends on the temperature of the elec- 
tronic reservoirs! 27 * 28 ^ The damping of the bosonic system 
at zero temperature, with rate Kp^ has also been taken 
into account by including the last Lindblad term in equa- 
tion ^. 

We are interested on the generation of qubit-qubit en- 
tanglement and on the transport properties in the sta- 
tionary state, p stat . This can be obtained from equa- 
tion ^ as p(t) = £p stat = such that the Liouvil- 



lian C has a zero eigenvalue with right eigenvector de- 



noted as |0)) 



The corresponding left eigenvec- 



tor is |0)) such that the probability conservation reads 
<<0|0>> = Tr[lp stat ] = 1. Using this language, the av- 
erage of any operator A acting on the qubits-resonator 
system reads (A) = Tr[Ap stat ] = «6|i|0}) = ((A)). 

The set of equations for the elements of the density 
matrix p nm (t), in the basis given by the direct prod- 
uct of the electronic states and the oscillator Fock states 
|cci,a2) ® \n) (with n = 0,1,2,...), is solved numeri- 
cally by truncating up to a maximum number of pho- 
ton states n — N max . We take the order of magnitude 
of the parameters from the recent experiments report- 
ing circuit-QED devices with semiconducting QDsP^^ 
Even though we focus here on this moderate coupling 
regime g/oJ r ~ 10~ 2 , we note in passing that our numer- 
ical scheme allows in principle to include stronger cou- 
plings, such as the ones already achieved i n cir cuit-QED 
architectures with superconducting qubitsP^Sl 



In practice, the shot noise at zero-frequency is calcu- 
lated in terms of the inverse of the part of the Liouvil- 
lian that is non-singular at zero-frequency (or pseudo- 
inverse), R = QC' X Q (with Q = 1 - |0))((0|), see 
e.g. Refs. HUH32J The diagonal part of the noise reads 
S«(0) = 2(((Ii)) - 2((liRIi))), with i = 1,2, whereas 
the off-diagonal noise cross-correlations read 6*12(0) = 
521(0) = -2 (((1 1 RI 2 )) + ((I2RI1))). Note that any fi- 
nite off-diagonal noise in this setup indicates correlations 
between distant currents across each DQD. 

In what follows we present our noise results in the form 
of Fano factors, defined as Fij = iSy(0)/(2e-\/l,-/j), which 
quantifies deviations from the Poissonian noise originated 
by uncorrelated carriers. In particular, super-Poissonian 
noise (F > 1) is related to a bunching behavior, of the 
carriers whereas sub- Poissonian noise (F < 1) signals 
ant i-bunching S3 



A. A single DQD coupled to the resonator 



III. STATIONARY TRANSPORT PROPERTIES: 

CURRENT, SHOT NOISE AND CURRENT 

CORRELATIONS 



We expect that the indirect, non-local two-qubit inter- 
action induced by the coupling to a common resonator 
mode can be revealed in transport through either DQD. 
As previously mentioned, we restrict ourselves to the 
Coulomb blockade regime in the infinite bias voltage 
limit. 

In this case of unidirectional transport, the total cur- 
rent passing through the DQDi is described by the 
operator Zj = eT^idn,ipd Ri , and the corresponding 
steady-state expectation value reads Ii — ((0|Zj|0)) = 
Tr[lip stat ]. 

We also analyze the non-equilibrium quantum noise, 
resulting from the temporal fluctuations of the current, 
by means of the current-current correlation function 
(A2i(r), AJj(O)), with AIi(t) = I t (t) - (F). The Fourier 
transform of such correlation function defines the power 
spectral density of shot noise: 



Sij(u) 



dre^dAI^AIM}) (6) 



It has been shown that this finite-frequency power spec- 
tral density contains a great deal of information about 
internal dynamics of the system!^ Nevertheless, we here 
restrict the analysis to zero frequencies for simplicity. In 
particular we focus on the cross-correlations which, as 
we shall show, exhibit features related with the qubit- 
qubit effective interaction induced by the common cou- 
pling to the resonator. Additional interest in studying 
shot noise and cross-correlations reside in theoretical pro- 
posals which make use of current correlations to study 
and detect entanglement in mesoscopic systems! 



To set the stage for our study, we begin by analyz- 
ing the case of a single DQD coupled to the cavity. The 
physics here is that of inelastic transport through a two- 
level system, a problem which has received a lot of atten- 
tion in various contexts! 28 * 44 " ^ In the frame of circuit- 
QED with semiconducting qubits, the problem has been 
theoretically studied in Refs. 19 and 20 mainly with focus 
on lasing. 

In figure |2k) we show the current in the DQD as a 
function of its level detuning E\ (all the parameters are 
expressed in units of the resonator frequency ui r ). As 
expected, there is an elastic peak in the current around 
£\ = which corresponds to resonant tunneling across 
the DQD. Here, the electronic transport occurs by the 
tunnel coupling with the reservoirs, which we assumed 
to be the same for both leads Tl,i = r^,i. The height 
and width of the clastic peak is in agreement with the 
well known analytical expression for the current through 
a DQDp2E3For finite detuning (i.e., with the electronic 
levels of the DQD far from resonance) the current is sup- 
pressed except at values of e\ corresponding to a res- 
onance condition at which the frequency of the qubit 
fii = \Je\ + 4i 2 equals the frequency of the resonator uj r . 
This feature corresponds to inelastic processes in which 
the tunneling of an electron between the left and right 
dots of the qubit excites the state of the resonator. This 
behavior is in qualitative agreement with the theoretical 
results of Jin et al. p3 who studied lasing in a DQD-based 
circuit-QED system (the main idea being that transport 
of electrons through the artificial two-level system can 
lead to a population inversion and induce a lasing state 
in the microwave resonator). Although we are not in- 
terested on analyzing the specific lasing conditions, the 
underlying mechanism giving rise to the inelastic peak of 
the current is the same.— Finally, we mention in pass- 
ing that the physics is also similar to one of spontaneous 
emission of a DQD coupled to a bath of phonons, demon- 
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FIG. 2: (Color online) Results for a) stationary current and 
b) Fano factor for a single-qubit coupled to a transmission 
line resonator. The solid line corresponds to equal tunneling 
rates to the reservoirs (r_L,i = r_R,i = 10 -3 ) and the dashed 
line to asymmetric rates (Fl.i = 0.01, Tu,! = 10~ 3 ). Note 
that the inelastic peaks appear at values of si corresponding 
to the resonance condition Q\ = nhus r . Inset: zoom of the 
current peak at fii = 2hw r for symmetric rates. The rest of 
the parameters (in units of uj r ) are: ti = 0.2, g\ = 0.008, 
k = 5x10" 4 . 



strated experimentally in Ref. [55] as well as the physics 
of on-chip noise detection using two-level systems! 45 * 52 * 56 ^ 

For large enough electron-boson coupling g%, addi- 
tional resonances at f2i rj nHu) r appear. An example 
for n — 2 is shown in the inset in figure [2k). 

Note that the qubit-photon resonances are not present 
for the region e\ < since the photon mode is in its 
ground state at zero temperature and, hence, photon ab- 
sorption inducing tunneling of the charge is not possible. 
As we shall discuss in the next subsection, this is no 
longer true in the presence of a second qubit that can 
induce excitations in the microwave resonator. 

The corresponding Fano factor Fu, shown in fig- 
ure [2b), exhibits a dip around £\ = 0. There, interdot 
tunneling delocalizes the charge which, combined with 
the strong Coulomb blockade, reduc es the noise and gives 
sub-Possonian Fano factor, Fu < l) 46 l 57 l As the level de- 
tuning £\ increases, the charge becomes localized, say in 
the left dot for E\ > 0, and hence Poissonian noise from a 
single barrier (the one parametrized by Tl) is obtained. 
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FIG. 3: (Color online) Steady state current in the first DQD 
as a function of its level position £i for two different con- 
figurations of the second qubit: £2 = —1 and £2 = 1- The 
two examples have been vertically shifted (with an offset of 
0.4) for the sake of clarity. Parameters (in units of uj r ): 
9i = 32 = 9 = 0.008, ti = t 2 = t = 0.2, T = 10 -3 and 
«= 10" 3 . 



This is so until the resonance conditions fii = nHuj r 
are reached, where the noise is reduced again yielding 
Fn < 1. This sub-Poissonian value at resonance with 
the photon cavity reveals that the charge is transferred 
across the DQD with the simultaneous excitation of the 
resonator. The same kind of result is obtained for emis- 
sion into a full bath of bosons. 4? 

Note also the small resonant feature in the region e\ < 
0. Even though in this configuration the extra charge is 
mainly localized in the left dot, there is a small probabil- 
ity of populating the right dot (and subsequently tunnel 
out from the right barrier). From the point of view of 
the qubit, this means that there is a small probability of 
populating the excited state and hence to emit photons. 
This can be easily seen if we write the qubit-photon in- 
teraction in the qubit eigenbasis |e) = cos ||i)+sin ||i?) 
and \g) = sin||£) - cos||i?), with 9 = arctan(^L) be- 
ing the angle that characterizes mixing in the charge 



1 AT* ) (at 



with 



subspace: H e _ res = gi(cos 9t z - 
T z = |e)(e| - \g)(g\ and t x = \e)(g\ + \g)(e\. We have 
checked that photon emission at £\ « — 1 is small but 
finite (the photon occupation has a resonance around 
this detuning and increases from zero to (n) « 1CP 3 , 
not shown), as a result of |e) — > \g) relaxation processes. 
Dynamically, these rare events, where the qubit is ex- 
cited for negative detuning such that photon emission is 
possible, contribute to the noise which shows a feature 
at f^i — u> r with £1 < 0. On the contrary, they do not 
change significantly the average current, demonstrating 
the superior sensitivity that noise has. 

The effect on the transport properties of asymmetric 
tunneling rates is also shown in figure [2l where we consid- 
ered that r^i > r.R,i- The current, figure [2k), exhibits 
the same qualitative behavior than the case with equal 
rates, with an elastic peak around e = and satellite 



peaks at the resonances qubit-resonator. On the con- 
trary, the Fano factor changes completely for asymmet- 
ric rates, figure^). In this case, i*n presents a double 
peak structure in the region e = 0, with the maximum of 
the peaks corresponding to super-Poissonian noise. This 
well-known effect can be understood from the analytical 
expression of the Fano factor^ and, in particular, is orig- 
inated from the smaller coupling to the drain reservoir 
which, ultimately, makes the Coulomb interaction more 
effective and gives rise to bunching in transport with 
Fn > 1. The same kind of bunching behavior is observed 
for the Fano factor at the qubit-photon resonances. It is 
interesting to compare this F\\ > 1 at the one photon 
resonance with the result for a full b osoni c bath which 
always results in sub-Poissonian noise JiSlSSl Hence super- 
Poissonian noise results from the qubit-photon coherent 
interaction. This result is also along the lines of Ref. 
[59l where the authors demonstrate that the bunching 
effect cannot be obtained from a picture without qubit 
coherences. In the context of lasing, this sort of super- 
Poissonian noise has been related to squeezing of the res- 
onator stateP3 



B. Two DQDs coupled to the transmission line 
resonator 

We turn now to our original model in which two DQDs 
are coupled to the same photon mode of the microwave 
resonator, but uncoupled to each other. For simplic- 
ity, we consider first the same intra-dot tunnel couplings 
ti = t and equal electron-photon coupling gi = g. It is 
assumed that the tunneling rates to left and right reser- 
voirs are equal and also equivalent for both DQDs i.e., 
^L,i = Fn,i = T, unless otherwise stated. As in the case 
for a single DQD, all the parameters are given in terms 



of the bare frequency of the microwave resonator ui r . 

Results for the mean value of the stationary current 
passing through the first DQD, ii, as function of its level 
detuning e\, while keeping the second DQD in a fixed 
level structure, are presented in figure [3] Similarly to 
the case for a single DQD, there is an elastic peak in 
I\ around £\ = 0. A second, inelastic peak appears in 
the region where this qubit enters in resonance with the 
photon, fii s=s hhj r revealing that this effect is entirely 
due to the coupling of this DQD with the resonator and 
thus will appear irrespective of the presence of the second 
DQD. We refer to this feature as the one-qubit (lqb) peak. 

Interestingly, an additional peak in I\ is observed in 
the emission part e± > for the case £2 = —1, and in 
the absorption part e\ < with e<i = 1. The third peak 
arises when both DQDs are brought in resonance with 
each other, fii = O2, with opposite detuning, e\ = —£2, 
but slightly out of resonance with the photon mode, 
fii = O2 7^ hio r . It is a result of an indirect qubit-qubit 
interaction induced by the common coupling to the mi- 
crowave resonator and therefore we refer to it as the two- 
qubits (2qb) peak. The fact that this resonance appears 
at an energy larger than the bare frequency u> r reveals 
that the effective interaction is obtained via virtual pho- 
tons: when both qubits are in resonance, the excitation 
in one of the DQDs is transferred to the other by virtually 
becoming a photon in the microwave resonator. Similar 
physics has been demonstrated experimentally in circuit- 
QED systems with superconducting qubits, see Ref. ItHl 

In order to have a better understanding of the induced 
qubit-qubit interaction, we derive an effective Hamil- 
tonian for the regime where the 2qb-features appear, 
fli — co r > g. Starting from the Hamiltonian Hq = 
H e i + Hint + H res + H e _ res , and restricting ourselves 
to states with n = 0, 1, we obtain to second order in the 
qubit-photon coupling: 



H eff = 



E ( 9 £ j°i + heff°i ) + J^Wz ~ E J xz,ij<rl<4 



E 



fi 2 
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1 



i l<i — LO r l li -\- k-V 






(7) 



The effective Hamiltonian of equation fff\ explicitly 
shows that the interaction of the qubits with a com- 
mon photon mode translates into a shift of their fre- 
quencies, through the renormalized tunneling ampli- 
tude U -> ti.pjt = U 1 



eff 



ML 



as 



well as two types of qubit-qubit interaction. The first 
one is Ising like with effective exchange constant J z = 



E 9192 
% n? 



2f 



( 1 

\Qi-Ulr. 



Qi+Ulr 



whereas the sec- 



ond one is an XZ exchange interaction with a coupling 



strength J xz 



9i9j£. 






where 



j \ J 3 r 

the presence of an electron in one qubit induces a tran- 
sition in the other. 



Neglecting the terms of order l/ui r and l/(fij+w r ) the 
effective Hamiltonian reads (the last, constant, term in 





1 ' 1 




1 


0.02- 


Eff. model (7) 

Eff. model (8) 

Full Master Eq. (5) 


1 


0.01 - 






■ 


0.00- 




V_^— 



0.94 



0.96 



0.98 

e./ox, 



1.00 



1.02 



0.08 



0.06 



0.04 









1 








e 2 =-0.98 










- - -E 2 =-1 






i 




E 2 =-1.1 












- 



0.98 1.00 



1.02 



1.04 
e,/(B, 



1.06 1.08 



1.10 



FIG. 4: (Color online) Comparison of the results for the 
steady state current in the first DQD as a function of its 
level position ei, obtained with the full Master Equation ([5| 
and with the models described by the effective Hamiltonians 
(JtJ) and J8j. Parameters: e 2 = -1, ffi = Q2 = g — 0.008, 
h = t 2 = t = 0.2, T = f0~ 3 and k = 10" 3 . 
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show dispersive shifts due to the off-resonant interaction 
with the microwave resonator photons.^ 

The interaction terms in equation Q capture quite 
well the 2qb transport resonance. This is explicitly shown 
in figure [4] where we plot a comparison of the current cal- 
culated with an effective master equation obtained from 
the models (JTl) and p| against the one obtained with 
the full master equation given by (kM), around the res- 
onance ill ~ ^2- It can be noticed that the effective 
models reproduce the width and height of the 2qb peak. 
Of course, outside this resonance the effective model fails 
and cannot describe transport in the full regime of level 
detunings. 

Once we have shown that the 2qb feature comes in- 
deed from a resonator-induced interaction between both 
charge qubits, we describe how the non-local character 
of this interaction can be easily explored. 

This is explicitly demonstrated in figure [5] where we 
show results for I x as a function of e\ around the two- 
qubit resonance condition and for different values of £2. 



FIG. 5: (Color online) Current on the first DQD as a function 
of £1 around the qubit-qubit resonance, for different values of 
£2. Rest of the parameters as those used in Fig. [3] 



The 2qb-peak in the current through one qubit clearly 
moves as one varies the level position in the other, while 
the lqb resonance remains unaltered (not shown) upon 
changing £2- We can also note that as the difference 
|Oj — uj r \ increases, the effective couplings given by equa- 
tion (pi) decrease and therefore the induced qubit-qubit 
interaction is turned off. Experiments along these lines 
have been recently reported for transport through single 
carbon-nanotubc quantum dots, where non-local control 
mediated by a photon cavity (in the classical limit) has 
been demonstrated.^ Thus we expect that an experi- 
mental test of our prediction in figure [5] is within reach. 

An even more interesting possibility is to non-locally 
manipulate the qubit-qubit induced interaction by tun- 
ing the dissipative coupling of one of the qubits with its 
fermionic reservoirs. For example, a strong coupling to 
the right reservoir in, say, qubit 2 induces a transport 
version of the quantum Zcno effect which tends to freeze 
the dynamics of the second qubit by effectively localizing 
the charge in the left dot of the DQD2, with (a^) —> 1. 
We demonstrate this effect in figure [6] where the current 
through the first DQD as a function of T^.2 is shown 
for the two-qubit resonance condition e\ — —£2, with 
r,3,i = ^l.2 = T. It is observed there that for fixed 
qubits parameters, the current through DQDI is strongly 
reduced by increasing merely the rate T^ of the second 
DQD. We can reinforce the interpretation of this results 
by recalling the effective qubit-qubit interaction: for very 
large Ffl j2 one can replace the operators of the second 
qubit by the corresponding mean value; then the effec- 
tive coupling constants for the first qubit are also frozen 
and results in a smaller effective coupling. 

To check more critically the presence of non-local cor- 
relations mediated by the cavity, we study how nontriv- 
ial noise correlations develop. The Fano factor for the 
DQDI shows the same qualitative behavior exhibited in 
the single-qubit case (for symmetric rates with the reser- 
voirs), with sub-Poissonian regions around all resonances 
of the problem, see figure [7j 
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FIG. 6: (Color online) Current on the first DQD as a function 
of Tft t 2 for the resonance s\ = —£2 = 1. Rest of the param- 
eters (in units of u r ): t = 0.2, V L ,i = 1^,1 = T L>2 = 1CT 3 , 
K= 10" 3 . 
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FIG. 8: (Color online) Correlators for the current passing 
through both qubits F21, for £2 = —1 and £2 = 1. Same 
parameters as in Fig[3] 
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FIG. 7: (Color online) Fano factor for the first DQD as a 
function of the level position ei for two different configura- 
tions of the second qubit: £2 = — 1 and £2 = 1. The results 
are vertically shifted (offset of 0.65) for the sake of clarity. 
Rest of the parameters as in Fig. [3] 



FIG. 9: (Color online) Colormap of the cross-correlations F12 
around the two-qubit resonance as a function of the level de- 
tuning £1 and the coupling parameter with the resonator, g. 
Rest of the parameters (in units of w r ): £2 = — 1, t = 0.2, 



r = 10" 



10" 



More crucially, the cross-correlations between sepa- 
rate currents through both DQDs, F12, develop sharp 
resonances at the qubit-qubit resonance f2i = H2 (fig- 
ure [8]). Apart from these clear resonances, other small 
features signal finite microwave resonator occupations 
which lead to non-zero correlations. As the coupling with 
the resonator increases, such features, and more gener- 
ally the overall behavior as a function of level detuning, 
can become rather intricate. Figure [9] shows the cross- 
correlations for increasing g in the region around the 2qb 
resonance. This figure reveals that the peak emerged 
around this resonance splits as the qubits-resonator cou- 
pling becomes larger. At the same time, the resonances 
become broader such that the function Fi 2 {£i,g) devel- 
ops a two-lobe structure. As we shall show in the next 
Section, this characteristic structure signals the forma- 
tion of Bell states between both qubits and hence the 
development of non-local entanglement. 



IV. QUBIT-QUBIT ENTANGLEMENT 

So far we have demonstrated that transport exhibits 
signatures of the induced interaction between the DQDs 
due to the common coupling to photons in the microwave 
resonator. Here, we go a step further an explicitly demon- 
strate that this common coupling can generate entangle- 
ment. In particular, we show that qubit-qubit entan- 
glement under nonequilibrium conditions can be gener- 
ated by virtual photons. For quantifying the nonequilib- 
rium entanglement we make use of the Concurrence,^ 
a measure of entanglement which arises from the entan- 
glement of formation and that is calculated by means of 
the density matrix of the system in the computational 
basis. We calculate the Concurrence of the steady state 
Pp atat , which corresponds to the projection of the sta- 
tionary density matrix onto the two-qubits subspace with 
a proper normalizationpZl and trace out the states of the 
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FIG. 10: (Color online) Concurrence for £2 = —1 and £2 = 1- 
Rest of the parameters as in Fig[3J 




0.98 0.985 0.99 0.995 1 1005 1.01 1.015 1.02 



FIG. 11: (Color online) Colormap of the Concurrence around 
the two-qubit resonance as a function of the level detuning 
£i and the coupling parameter with the resonator, g. Rest of 
the parameters (in units of uj r ): £2 = — 1, t = 0.2, T = 10 -3 , 
K= 10~ 3 . 



bosonic mode. 

Numerical results for the Concurrence, C, considering 
the same interdot tunneling amplitude in both qubits, 
ti = t, and symmetric electron-photon coupling gi = g 
are shown in figure [lO] for two different level detunings in 
the second qubit. For the typical value of the coupling 
g = 0.008 used here, C shows sharp features in the 2qb 
resonance, £1 = — £2. 

In figure [IT] we show the detail of Concurrence in the 
region of the 2qb resonance, as a function of £1 and 
the coupling strength to the microwave resonator g, for 
£ 2 = — 1. Here we find that, in the same way as the 
cross-correlators F±2 (figure [9]), the peak exhibited by 
the Concurrence around resonance splits and develops a 
two-lobe structure as the coupling g becomes larger. The 
similarity between these two quantities shows that cur- 
rent cross-correlations in the above configuration consti- 
tute an indicator of non-local qubit-qubit entanglement. 

The previous interpretation is supported by an analy- 




FIG. 12: (Color online) Stationary occupation probability of 
the Bell state \4> + ) (left) and \<j>~) (right) around the two- 
qubit resonance. Same axes and parameters as in Fig. [TTJ 



sis of the steady state populations of the system. If the 
analysis if done in terms of the populations in the local 
basis (e.g. |ai,cv2), with a = L,R) the double-peaked 
structure of Figs. [TTJ and [9] is hard to explain, since all 
local populations exhibit just a single peak around reso- 
nance. However, considering the stationary populations 
in the Bell basis of maximally entangled states, a different 
picture arises. Figure fl2] shows the population of the Bell 
states |<I> ± ),I^ which written in terms of the occupation 
of the L/R dots of each DQD read: 



14^ = 



1 



\R l ,R 2 )±\L 1 ,L 2 )) 



(10) 



The occupation probability of these two states show a 
double peak structure as g becomes larger. Importantly, 
these peaks occur asymmetrically such that each Bell 
state has maximum occupation on either side of the res- 
onance. The two-lobe structure of both the Concurrence 
and the cross-correlations thus correspond to the two 
maxima of the 1$*) Bell state populations. The remain- 
ing Bell states 



1^ 



1 



V2 



{\R 1 ,L 2 )±\L 1 ,R 2 )) 



(11) 



just show single peaks centered on resonance and pre- 
sumably do not greatly contribute to the overall form of 
the current cross-correlations. 



V. ASYMMETRIC COUPLING TO THE 
MICROWAVE RESONATOR, gi / g 2 . 

Finally, we also explore the effect of asymmetric val- 
ues of the electron-photon coupling strengths for each 
qubit, g\ ^ g 2 . Experimentally, this asymmetry can be 
achieved by changing both the capacitive coupling of each 
DQD to the microwave resonator Cf as well as the capac- 
itance of each DQD to ground C® , as the couplings scale 

; rP Our motivation here is to explore the 



as gi 



Cf+Ci 



possibility of detecting the interaction-induced shifts di- 
rectly in transport. Further motivation comes from Ref. 
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FIG. 13: (Color online) Colormap of the current around the 
qubit-qubit resonance condition as a function of the level de- 
tuning E\ and the coupling parameter with the resonator for 
the second qubit, <?2, for fixed g\ — 0.008. The dashed line 
indicate the values of e\ corresponding to the renormalized 
frequency of the first qubit Cli, e ff = ^2,e// predicted by the 
effective Hamiltonian pj. Rest of the parameters: £2 — — 1, 
£ = 0.2, T = 10" 3 , k= 10~ 3 . 



I66[ which theoretically proposed the use of inhomoge- 
neous coupling between two-level systems and a single 
quantized mode to generate and control multipartite en- 
tangled states 

The current as function of s± and gi is shown in fig- 

For increas- 



13 for the region around E\ 



-£2- 



ing §2, the position of the resonance is shifted with re- 
spect to the initial value for g\ — g 2 - As expected, 
this can be understood by means of the renormaliza- 
tion of the intra-dot tunneling coupling £j, e // in the ef- 
fective Hamiltonian of equation (fSJ) . This renormaliza- 
tion leads in turn to a change in the frequency of the 

qubits as fii, e // = J £J ~^^ieff Therefore, the current 
shows a dispersive shift at values of E\ accordingly to 
the new, effective qubit-qubit resonance condition given 
by £li <e ff — Q,2,eff- The dispersive shift obtained with 
the full numerics agrees with the one given by the ef- 
fective Hamiltonian, represented by the dashed line in 
figure [l3j Measurements along these lines would con- 
stitute further proof of resonator-induced interaction be- 
tween qubits. The same dispersive shift is also observed 



in the shot noise cross-correlations, figure 14 1), where 



again, the 2qb resonance in F12 splits for large enough 
coupling. 

Finally we present the same analysis for the Concur- 
rence in figure n4p). Apart from the shift, we can no- 
tice that, in general, the Concurrence has larger values 
in comparison to the case with g\ — g2, indicating that 
the asymmetry between the coupling parameters of each 
qubit with the bosonic mode makes the qubit-qubit en- 
tanglement to be more robust. 




FIG. 14: (Color online) Colormap of a) cross-correlations 
-F12 and b) Concurrence as a function of e% and (72, for 
fixed g\ = 0.008, and around the two-qubit resonance. The 
dashed line indicate the renormalized two-qubit resonance 
^i,e// = ^2,e//- Rest of the parameters as in Fig. [13] 



VI. CONCLUSION 

We studied theoretically photon-mediated transport 
and the generation of steady state correlations between 
two open charge qubits defined in spatially-separated 
double quantum dots which are coupled to a common 
transmission line resonator. Our results demonstrate 
that the qubits are entangled due to the indirect coupling 
induced by photons in the microwave resonator. Consid- 
ering that each qubit is open to electronic reservoirs, we 
have analyzed their transport properties and found that 
they reveal the qubit-qubit interaction. In particular, we 
calculated the zero-frequency shot noise and the current 
cross-correlations as a function of the level detuning of 
one of the qubits, and observed the presence of different 
resonant features in the regions where the qubit enters in 
resonance with the photon as well as with the other qubit. 
In the examples we studied here, the quantum correla- 
tions involved in the transport of charge and which are 
responsible of the signal in the cross-correlations, yield 
in a finite value for the Concurrence when the qubits 
interact due to off-resonant photons. Therefore, we pro- 
pose that measurements of current correlations could be 
used as a possible method for detecting entanglement 
and, in general, qubit-qubit interactions mediated by the 
microwave resonator. This proposal is motivated also in 
the context of recent experimental achievements demon- 
strating the coupling of semiconductor QDs to microwave 
resonators.™^ 

The model presented here constitute a step further in 
the study of this kind of hybrid systems, which can be 
relatively easily extended to several qubits. In general, 
this system let us to explore the interplay between co- 
herent interactions, entanglement and the effect of dis- 
sipation and noise. Moreover, our model can also be 
applied to charge qubits defined by Cooper-pair boxes or 
to systems in which the quantum dots are coupled to a 
nanoelectromechanical resonator. 

Note added: While finishing this manuscript, two the- 
oretical works analyzing a similar setup to ours have ap- 
peared. The first work 68 focuses on the effect that non- 
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local interaction between two DQDs resonantly coupled 
to the oscillator has on finite bias voltage transport prop- 
erties (which are prone to finite temperature effects in 
the electronic reservoirs). In contrast, we here focus on 
a different operating regime where the non-local inter- 
action is induced off-resonance and transport occurs at 
very large voltages. In this large-voltage regime, the re- 
sults are essentially independent on the electronic reser- 
voir temperature and are valid at arbitrary couplings to 
the reservoirs. This large voltage regime is also analyzed 
in the second workpSl where some overlapping results 
about photon-mediated transport and finite shot noise 
cross-correlations have been reported. 
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